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We confront the recent improved measurements of the Dj^) (semi) leptonic decays with Lattice QCD (LQCD) with Uf = 
3 flavours and QCD spectral sum rules (QSSR) predictions. D -^ /ii^^ leptonic width data compared with theoretical 
determinations of fo, leads to the value of the CKM mixing angle : \Vcd\ = 0.230 (lO)oxp (9)th- Measured ratio of the 
D semi-leptonic widths combined with LQCD and QSSR predictions leads to the average |Vcd|/|V'cs| = 0.2175(88), and 
then to |Vcs| — 1.068(47). We consider the previous determinations as improvements of the existing estimates. Using 
the average data of the Ds -^ ^v^ (resp. Ds — > tVt) branching ratios, one obtains : |Vcs|/^ = (259 ± 12) MeV (resp) 
l^csl/i) ~ (274 ± 13) MeV, which can be compared with the average of the Standard Model (SM) values from LQCD and 
QSSR fos ~ (240 ± 7) MeV. If one uses the present determination of IVcal, there is an agreement with the SM prediction 
within la. If instead, we impose the unitarity constraint \Vcs\ < 1, or assume (as frequently done) |T4s| = \Vud\, we would 
obtain a deviation from the SM expectations ranging from 1.5 to 3 a, therefore, signaling some New Physics beyond the 
SM. 



1. Introduction 



Good experimental and theoretical controls of the lep- 
tonic decay constants fp (analogue to /^) of charmed 
D(^s) Bind beautiful i3(s) (pseudo)scalar mesons are of 
prime importance for understanding the dynamics of 
the heavy light quark systems (overlap of the wave- 
functions, heavy quark mass-behaviour,...) in QCD. 
Many efforts have been devoted to this study both the- 
oretically and experimentally Q The leptonic partial 



width of the -Dj" is normalized as: 



T{D+ ^l+^i)^^\V,,\'flmfMn, 



1 



ML 



(1) 



where q = d,s; Mn^ and mi are respectively the 
masses of the £>+ meson and of the /+ charged lep- 
ton. Gp is the Fermi constant and |Vcg| is the CKM 
matrix element controlling the weak coupling of the 
c and q quarks, fn^ is the leptonic decay constant 
normalized like /^ = (130.4 ± 0.2) MeV as: 



{rrig 



'\q-/5c\D,) = fo^Ml 



(2) 



The previous leptonic process is helicity suppressed 
(~ mf) and the last factor in parenthesis is the familiar 
phase space factor. While the electron mode is tiny 
due to strong helicity suppression, the measurement 
of the T mode involves the detection of additionnal 
neutrinos, such that, the muon mode is experimentally 
the cleanest and most accessible one. Therefore, we 
shall mainly consider the D^ —> ijl^v^ data which (in 
principle) provide the most precise measurements, and 
show for a comparison the one from Ds — ^ tVt- decay. 

2. Determination of \Vcd\ 

For determining fu from D — > /ij/^ decay, it has be- 
come usual to assume that | Vcd \ = | Kis | , which is an 
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[T]; For earlier reviews, see e.g. 



extra input in the analysis. In addition, the extrac- 
tion of I Vus I from Kl-^ is affected by the uncertainties 
on the value of the K-form factor, while r-decay data 
give a slightly different (though consistent) value [3]. 



Table 1 



fn 



and fDs/fo from LQCD with nj — 3 flavours and QSSR. 



Method 



fo 



fnjfn 



lqcdH] 

LQCDi] 

QSSR 

Full QCD [213171 
Full QCD [8] 
Analytic Cont.jg] 



201 ± 17.3 
208 ±4 



203 ± 20 
195 ± 20 



249 ± 16.3 
241 ± 3 



235 ± 24 



1.24 ±0.07 
1.162 ±0.009 



1.15 ±0.04 



1.16 ±0.03 



Average 



202 ±8.3 



241.7 ±9.7 1.178 ±0.022 



The good agreeement between the measured and the- 
oretical values of fu [IIIOIII] . allows us to perform 
(with a good confidence) an opposite, though natural, 
procedure, by relying on the theoretical value of fo 
for extracting more accurately \Vcd\- By combining 
the most recent measured D -^ /if^ branching ratio 
from CLEO pUITT] : 



Br{D+ 



/u+i/^ 



) = (3.82 ±0.32 ±0.09) x 10" 



(3) 



with the averaged value from LQCD and QSSR a la 
SVZ 51] given in Table [TH, and using the D+ lifetime 



•^We have only considered the ratio obtained in [9] as the decay 
constants obtained there are systematically lower than other 
predictions, which may question the systematics of the Analytic 
Continuation method ^, that can cancel in the ratio. We have 
only quoted the most recent 3-loop predictions relevant for our 
discussions. More complete references are given in the reviews 
|2I3| . We have not taken the weighted average which should be 
dominated by the most accurate LQCD determination where the 
validity of the errors are still under discussions |13| . Instead, 
we have taken a nai've average of the central value and have 
averaged the errors quadratically. 



of (1040 ± 7) ps [H], we deduce: 

\Vcd\ = 0.230 ± O.OlOexp ± 0.009th , (4) 

to be compared with : 

|Ks| = 0.2165 ±0.0026cxp± 0.0005th (5) 

from T-decay [4]. 

3. Determination of \Vcs\ 

Semileptonic decays of the hght K and heavy B 
mesons have given most of the CKM mixing param- 
eters. In the case of the process D^ -^ P~l^vi, 
{P^ = -K^ and Ps = K^) and [I = e, /i), where one 
can neglect the ratio {mi/rricY, the differential rate 
can be expressed as: 



dT 



D«^p-l+iyi 



^pl\Vc,\vr'^{i-')\-' 



rD^P /„2\\2 



(6) 



where q^ is the invariant mass squared of the Iv system; 
f^^^{q^) is the form factor and Pq is the momentum 
of the light meson. CLEO data lead to the result [15,: 



|/^^^(0)||v;.| 



= 0.188 (8) (2) , 



(7) 



where the errors are respectively statistical and sys- 
tematical. In order to use this data, we update the 
value of the ratio of form factors from QSSR [TB] by 
using the recent value ?fis(2 GeV) = (96.3 ± 17.5) MeV 
[T7] instead of the one ms{2 GeV) = 115.7 MeV used 
in the original paper. We have also used the value of 
rric in [2l3j Q Then, we obtain: 



ro ^ l/j^^^ (0)1/1/^^^(0)1 = 1.11 ± 0.07 , 



(8) 



which has almost the same value as the previous one. 
The error has been multiplied by V2 in order to take 
into account unknown higher order terms in the QCD 
series. 

Table 2 

Theoretical predictions of form factors. 



Methods 



l/f^^-'WI/l/f^'TOI 



LQCD [U 
LCSRlTSl 
QSSR[16] (updated) 



1.149 ±0.040 ± 0.119 
1.19 ±0.06 
1.11 ±0.07 



Average 



1.155 ±0.043 



which we compare in Table [2] with some other deter- 
minations, from which we deduce the average: 



TO ^ 1/^^^(0)1/1/^- 



'^(0)1 = 1.155 ±0.043 



(9) 



Some remarks are in order here: 

• The determinations of the ratio compared to the ab- 
solute values are always more accurate in the QSSR 



*There is an unfortunate misprint in |16| . Instead of 0.007 for 
the error in Eq. (24), one should read 0.07. 



calculations due to the cancellations of different sys- 
tematics. The sensitivity to the higher state contri- 
butions (continuum threshold-dependence), to PT ra- 
diative corrections (which tend to cancel in the ratio), 
to the exact value of the charm quark mass (absent 
to leading order) ,... are less pronounced in the direct 
evaluation of the ratio than in each individual form 
factor. The dependence on the sum rule variable of 
the prediction is also weaker as shown in the figure of 

• This value does not come from a numerical fit but 
from a semi-analytic expression where all different 
sources of corrections can be easily controlled. 

• The sum rule used here has successfully predicted 
(though with a modest accuracy) the absolute value of 
the form factor /f ^''(O) = (0.23±0.02) for B -^ n^iv^, 
and of other semileptonic decays [19] and the decay 
rate B^K*-f [20]. 

• As a good illustration of this property, the LQCD 
recent most precise prediction of the ratio [6] : 



/d,//i3- 1.162 ±0.009 



(10) 



agrees within a digit with the central value from QSSR 
prediction obtained one decade earlier |7|2|3j Q 



/d.//d- 1.15 ±0.04 



(11) 



which is a strong indication on the reliability of the 
QSSR approach despite its modest accuracy u, and a 
compatibility of the LQCD and QSSR results as often 
encountered in different channels. The same remark 
also holds for the absolute value of the decay constants. 
• It is also important to notice that the previous ratio 
in Eq. (|lip comes from a semi-analytic functional de- 
pendence, which is easy to control rather than from a 
numerical fit. 

Using the previous inputs and taking into account pre- 
vious remarks, we can deduce with a quite good accu- 
racy : 



Wr 



cd\ 



14.. 



0.2171 (95)e,p (81) 



jth 



(12) 



which is compared in Table [3] with some alternative 
determinations 15J. Using the previous value of \Vcd\ 
into Eq. p^ . one can deduce: 



|14s| = 1.059 (65)exp (58)th. 



(13) 



Considering as a final result, the (weighted) average of 
the three determinations in Table [31 we deduce: 



Wr 



cd\ 



w. 



= 0.2175 



114,1-1.068(47) 



(14) 



cs 



^An alternative estimate using the analytic continuation method 
leads to similar result but with a smaller error 9 as given in 
Tabled 

^As indicated in 7 2i3], the error in this ratio has been enlarged 
by \/2 for a crude estimate of the unknown m'^a^ corrections, 
which I wish will be available in the near future. However, after 
inspection of the analytic expression of the spectral function this 
crude estimate can give an overestimate of the true error. 



to be compared with the usual assumption: 



\Vc>_ 



\Vud\ =0.97377 ±0.00027. 



(15) 



We consider these results as improvements of earlier 
results obtained in [TB] . Further tests of the assump- 
tion \Vcs\ — I Kid I and some eventual deviations from 
the unitarity conditions can be reached in future im- 
proved measurements of _D-decays. These eventual de- 
viations may also reveal some New Physics beyond the 
SM expectations. 

Table 3 

IKdl/lVcsl&IVcsl using CLEO data + Theoretical mehods. 
Tile 1st (resp) 2nd errors are experimental (resp) theorctieal. 
Tiie 3rd error in CLEO is due to /°^^(0) from LQCD. 



Method \V^d\/\Vc, 



\V.A 



Comments 



LQCD[6] 0.2257(209)(20) 1.070(70)(10) 

CLEO[l5] 0.2138(100)(30)(110) 1.075(69)(47)(55) LQCD 
This work 0.2171(95)(81) 1.059(65)(58) QSSR 



Average 0.2175(88) 


1.068(47) 




Table 4 

Ds —> fJ:!^^ branehing 


ratios. 




Exp. 


S„ X lO-'^ 


B*,(%) 


CLEO-e[211 

BELLE[22l 

CLEO 23] 

BEATRICEI24] 

ALEPH[25l 

BABAR[26l 


5.94 ±0.66 ±0.31 
6.44 ±0.76 ±0.52 

6. 2 ±0.8 ±1.3 ±1.6 

8. 3 ±2. 3 ±0.6 ±2.1 
6.8 ±1.1 ±1.8 

6. 74 ±0.83 ±0.26 ±0.66 


3.6 ±0.9 
3.6 ±0.9 
3.6 ±0.9 
4.71 ±0.46 


Average 

(no rad. eorr.) 


6.13±0.57+> 
6.33 ±0.47 





+) Exclude <hn~^ mode normalizations. 



4. The value of fo^ from Dg 



^IVf, 



We shall use the previous value of the CKM angle 
\Vcs\ in Eq. (HH) for extracting the value of fn^ 
from Ds — > nVfi , which has been emphasized to be the 
cleanest experimental mode. We shall use the average 
of the different data in Table [4] and the Dg lifetime of 
(0.500 ± 0.004) ps [Mj. One can notice that the av- 
erage including or excluding the ^tt"*" modes is almost 
the same. However, taking into account that the nor- 
malization to the </)7r"'" modes induce more systematic 
uncertainties [T], we shall only consider the 1st aver- 
age which does not use this normalization. Radiative 
corrections will decrease the branching ratio by about 
2%, which we shall include in the extraction of the 
decay constant. Then, we deduce: 



/^ - (242.5 ± 10.7v.. ± 11.5exp) MeV 



(16) 



^In the current literature, one often assumes \Vcs\ = IV'udl- 



The uncertainty is comparable with the one from the 
theoretical average given in Table (TJ and is larger than 
the one quoted in [llllj . The main reason is that the 
result of |llllj uses the assumption \Vcs\ = |14d| (see 
Eq. p^ . where \Vud\ has a tiny error. The central value 
of f'^^ is also smaller than in |llll| because the one 
of \Vcs\ obtained in Eq. (fT^ is larger than that of 
\Vud\- If instead, we only impose the unitarity con- 
straint \Vcs\ < 1, one can deduce the lower bound: 



WrJ < 1 



/^ > (259±12) MeV (90% CL) .(17) 



The usual assumption \Vc 
would imply: 



IK. 



IK, 



/," 



Ds 



= \Vud\ used in [nn] 
(266 ±12) MeV. (18) 



5. The value of fo^ from Ds -^ tv^ 

Though expected to be less accurately measured than 
the -Ds -^ fJ-i^fj, decay, it is informative to repeat the 
previous analysis for the Ds — *■ tVt data, which are 
given in Table [HI 



Table 5 



Ds 



■ tUt branehing ratios. 



Exp. 



Br X 10^ 



CLEO-eim 
CLEO-el271 
ALEPH[25l 
L3[28l 
OPAL [29] 



8.0 ± 1.3 ±0.4 
6.17 ±0.71 ±0.36 
5.8 ± 0.8 ± 1.8 
7.4 ± 2.8 ± 1.6 ± 1.1 
7.0 ±2.1 ±2.0 



Average 



6.57 ±0.63 



From the averaged branching ratio, one can deduce in 
MevB 



fk 



256.6 ± 11. 3v,,±12.4o 



IK, I in Eq. m 



> 273.5 ±13.2 
= 280.9 ±13.5, 



cxp 



iKsl < 1 (90% CL) 



cxp 



IK,,I = IK 



ud\ 



(19) 



where one can notice that the value of fo^ from r-data 
is about Ict higher than the one from //-data: 



fLJfD =1.06 ±0.07 



(20) 



6. Comparison with theoretical predictions 

One can compare the previous experimental values 
with the theoretical predictions reviewed in 1 2 ''3], 
which we have selected in Table [T] From this table, 
the average of the direct determinations of /d^ is: 



f^'l ^ (241.7 ±9.7) MeV , 



(21) 



while the one obtained from the average value of Jd 
multiplied by the ratio fosl Id is: 



y-ratio ^ ^238.0 ± 11.8) MeV 



(22) 



*As mentioned in [T], the t rate does not need to be radiatively 
corrected. 



from which we deduce the nai've average : 



j-th 

Id, 



(240 ± 7) MeV 



(23) 



where the errors have been added quadraticaUy. We 
consider this result as the final theoretical result to be 
compared with experiments. 

• /Ds from Ds — > fii^^- decay: one can notice that the 
value of /d, obtained in Eq. P^ by using the value 
of |V"cs| in Eq. p^ agrees quite well with the previous 
average of the SM predictions in Eq. ([^5]) . Discrepan- 
cies between the experimental and theoretical numbers 
start to be visible (> 1.5 a) when using the unitarity 
constraint for \Vcs\ < 1 [ see Eq. (|17l) ] and becomes 2 a 
when one uses the additional assumption | Vcs \ = \ Vud \ 
[ see Eq. (UHl) ] . 

• /-Ds from Ds — > tVt -decay: here, a comparison of the 
theoretical prediction shows larger discrepancies rang- 
ing from 1 CT in the case of \Vcs\ from Eq. (|14p to 3 cr 
in the case iKsl = |Kd|- 



• First, these data have been used to estimate the 
CKM mixing angles \Vcd\ and \Vcs\- The results are 
given in Eq. ^ and in Eqs. ^ to pi)) . 

• The resulting value of \Vcs\ have been used together 
with the Ds — > /ii/^ and Ds — > tVt data for extracting 
the leptonic decay constant /^^ , which behaves like 
1/lK.I [Eqs. IHD to dH]. 

• Comparing the experimental value of fo^ with the 
theoretical averaged LQCD and QSSR determinations 
in Eq. ((23)) . we conclude that a sharp detection of an 
eventual deviation of the SM predictions from these 
processes, needs either better determinations of the 
CKM mixing angles \Vcd\ and \Vcs\ than the one ob- 
tained in this paper, or a strict validity of the fre- 
quently used assumption \Vcs\ = \Vud\ or of the uni- 
tarity condition \Vcs\ < 1- In these cases, the data 
indicate deviations from the SM expectations from 1.5 
to 3 a^ therefore, signaling some New Physics beyond 
the SM. 



7. Interpretations: signals of New Physics ? 

• Some deviations from the SM expectations can be 
manifest if the |Vcs| satisfies the unitarity constraints. 
Therefore, more precise determinations of \Vcs\ and to 
a lesser extent of \Vcd\ are required for a sharp test of 
the SM predictions. Some eventual deviations of these 
CKM matrix elements from unitarity constraints can 
also signal some departures from the SM expectations 
independently of the fu^ values. 

• We have shown that the deviation of fo^ [Eqs. p6|) 
to p9|) ] from the SM expectations [Eq. ((23l) can be in- 
deed quite large (3cr in the r channel) but a sharp 
conclusion needs a better control of the CKM mixing 
matrix [Ksl- 

• We have also found in Eq. (P0|) that the value of fo^ 
from Ds -^ tVt decay can be larger than the one from 
Ds -^ fii^fj. decay. More precise measurements of fu^ 
in both channels can then provide an universality test 
of the Wfii'^ and Wrvr vertices for these processes. 
At present, t decay data show a slight violation ^■. 



Bw^r.jBw-.^., = 1.039 ±0.013 



(24) 



which can permit some new interactions (not neces- 
sary the same) at these vertices. 

• Using the fact that the value of the decay constant 
fu^ deviates from the SM predictions, some attempts 
to explain this deviation as due to non-standard effec- 
tive interactions or/and to new scalar particles or/and 
to leptoquarks have been discussed in the literature 
[50] , Analysis of different models beyond the SM are 
beyond the scope of this paper though planned to be 
done in a future work. 

8. Conclusions 

We have re-examined the consequences of the recent 
data on (semi) leptonic Du) decays by confronting 
them with the LQCD and QSSR theoretical calcula- 
tions. 
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